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Abstract 

The  introduction  of  high-speed  circuits  to  realize  an  arithmetic  function  /  as  a  piecewise 
linear  approximation  has  created  a  need  to  understand  how  the  number  of  segments 
depends  on  the  interval  a  <  x  <  b  and  the  desired  approximation  error  e.  For  the  case 
of  optimum  non-uniform  segments,  we  show  that  the  number  of  segments  is  given  as 
s(e)  (e  -»•  0),  where  c  =  \  £  vlT"!  x)\dx.  We  also  show  that,  if  the  segments  have 

the  same  width  (to  reduce  circuit  complexity),  then  the  number  of  segments  is  given  by 

«(e)  ~  -»•  0),  where  c  = 
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1.  Introduction 

A  numeric  function  generator  (NFG)  is  a  logic  circuit  that  realizes  an  arithmetic 
function  like  f(x)  =  sin  x  over  some  specified  interval  a  <  x  <  b.  We  consider  the  numeric 
function  generator  shown  in  Fig.  1.  This  architecture  realizes  a  given  function  as  a  set 
of  segments  or  pieces,  where  /  is  approximated  in  each  segment  by  the  linear  equation 
f(x)  ~  c.\X  +  Co-  The  values  of  C\  and  Cq  are  stored  in  the  Coefficients  Memory  shown  in 
Fig.  1  in  a  location  whose  address  is  specified  by  the  Segment  Index  Encoder.  In  designing 
this  circuit,  one  partitions  the  interval  a  <  x  <  b  into  segments,  where  the  coefficients 
Ci  and  Co  are  the  same  in  each  segment  and  approximate  the  function  to  within  some 
specified  error.  It  is  known  [9]  that  the  Segment  Index  Encoder  is  tractably  realized. 

This  circuit  is  said  to  realize  a  non-uniform  segmentation  because,  in  general,  the 
segments  have  different  widths.  Indeed,  we  will  choose  the  segment  widths  as  large  as  pos¬ 
sible  so  that  the  approximation  is  still  less  than  or  equal  to  within  the  given  approximation 
error.  In  this  way,  we  produce  a  segmentation  that  has  a  few  segments  as  possible. 
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Figure  1.  Architecture  of  a  Numerical  Function  Generator  Using  Piecewise  Linear  Ap¬ 
proximation  and  Non-Uniform  Segmentation. 


In  contrast,  we  also  consider  a  uniform  segmentation.  In  this  case,  all  segments 
have  equal  width.  If  the  segment  widths  are  2m,  then  the  segment  index  encoder  can  be 
removed  and  the  higher  order  n  —  m  bits  used  to  drive  the  address  of  the  Segment  Index 
Encoder,  where  n  is  the  number  of  bits  to  encode  x.  In  this  case,  the  circuit  is  likely  to 
be  smaller  and  faster. 

Lip  to  this  point,  we  have  not  had  an  analytical  tool  to  predict  the  size  of  the  Coefficients 
Memory  as  a  function  of  the  function  realized,  the  domain,  and  the  required  approximation 
error.  That  is,  our  understanding  of  its  size  has  only  been  through  experimental  results; 

1. e.  specific  implementations  [5,7-10].  We  derive  an  expression  for  the  number  of  segments 
(size  of  the  Coefficients  Memory)  for  both  the  non-uniform  and  the  uniform  case. 

2.  Non-Uniform  Approximation  With  Unrestricted  Slope 

Let  /  be  a  three  times  continuously  differentiable  function  defined  on  the  domain  [a,  b\. 
In  the  case  of  unrestricted  slope,  our  algorithm  proceeds  by  generating  a  segmentation 
{xo,a:i, . . .  ,a;n}  of  [a,  b]  with  the  property  that,  in  each  of  the  n  segments  [xj,Xj+i],i  = 
0, . . . ,  n  —  1,  the  chord  between  xt  and  Xi+\  produces  a  linear  approximation  ( c\x  +  Co)  to 
/  within  a  previously  specified  approximation  error  £  : 

\f(x)  -  cix  -  c0|  <  £,  xe[xi,xi+i\. 

In  the  example  of  the  sin( nx)  function,  experimental  results  show  that,  for  a  specified 
approximation  error  £,  the  segmentation  algorithm,  in  the  case  of  unrestricted  slope, 
determines  a  number  of  segments  s  that  is  proportional  to  1  j \fe.  We  now  show  that  this 
is  a  general  result  for  a  large  set  of  functions. 

Specifically,  we  give  an  asymptotic  approximation  for  the  number  of  segments  s(e) 
needed  to  approximate  a  given  function  f(x)  to  within  a  given  approximation  error  e. 
We  say  that  t(e)  is  an  asymptotic  approximation  to  s(e),  expressed  as  s(e)  t(e),  if 
linu^0+  s(e)/f(e)  =  1.  In  our  use  of  this,  s(e)  is  the  exact  number  of  segments,  while  t(e) 
is  an  approximation  to  s(e),  which  has  a  simple  form.  Intuitively,  we  expect  the  number 
of  segments  to  increase  as  the  approximation  error  £  decreases.  We  seek  to  determine 
this  relationship,  since  it  provides  insight  into  how  hardware  complexity  depends  on  the 
approximation  error. 

In  what  follows,  we  divide  the  domain  [a,  b]  into  two  sets  depending  of  the  value  of 
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f'\x).  Let 

A£  =  {x  G  [a,  b]  :  \f"(x)\  <  and  (1) 

Be  =  {x  e  [a,b]  :  \f”(x)\  >  \/e},  (2) 

where  e  >  0.  Note  that  1)  A£  is  a  closed  set;  2)  Be  is  an  open  set,  3)  A£f]B£  =  0;  and  4) 
A£\JB£  =  [a,  b\ . 

Theorem  1.  Consider  a  piecewise  linear  approximation  of  f  on  the  domain  [a,b\  that 
is  accurate  to  within  e,  using  a  piecewise  linear  segmentation.  Let  f  be  three  times  con¬ 
tinuously  differentiable  on  [ a,b\ .  Then,  s(e),  the  number  of  segments  in  an  optimum 
segmentation  of[a,b],  satisfies  the  following  asymptotic  approximation: 

0),  (3) 


where 

c-\fa  V\f"(x)\dx- 

Proof  See  Appendix. 

For  example,  if  we  take  f(x)  =  sin7ra:  on  [0,  |],  we  find 


s(e)  ~ 


71 


J^2  \J sin  7tx  dx 
4  \[e 


Using  MAPLE  to  numerically  evaluate  the  integral  in  (5)  yields 

,  ,  0.2995 

4S)~  — 


(4) 


(5) 


The  following  corollary  of  Theorem  1  relaxes  the  conditions  on  the  end  points  of  the 
domain. 


Corollary  1.  Let  f  be  three  times  continuously  differentiable  on  the  open  interval  ( a,b ) 
and  \JJ71  be  improperly  Riemann  integrable2  on  the  closed  interval  [a,  b] ,  with  integrable 
singularities  at  the  endpoints  a  or  b.  Then,  s(e),  the  number  of  segments  in  an  optimum 
segmentation  of  [a,b\,  satisfies  the  following  asymptotic  approximation: 

0),  (6) 

V£ 

where 

c=\Ja  \/| f"ix)\dx.  (7) 

2The  function  /  is  improperly  Riemann  integrable  on  [a,  b]  if  /  is  Riemann  integrable  on  every  open 
subinterval  (c,  d)  of  [a,  b],  and  the  function  /  becomes  unbounded  in  neighborhoods  of  a  or  b ,  and 
limc^0jd^fc  /cd  f{x)dx  exists. 
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The  significance  of  Corollary  1  is  that  we  can  obtain  an  asymptotic  approximation  to 
the  number  of  segments  even  when  the  function  has  singularities  at  the  endpoints  of  the 
interval  over  which  the  approximation  occurs. 

Table  1  shows  the  asymptotic  approximations  to  the  number  of  segments  for  14  func¬ 
tions.  These  values  were  derived  from  Theorem  1  and  Corollary  1,  and  are  shown  in  the 
third  column  labeled  Non-uniform.  The  two  functions,  \fx  and  the  entropy  function, 
— x  log2  x  —  (1  —  x)  log2(l  —  x ),  on  which  we  applied  Corollary  1,  are  among  the  functions 
with  the  largest  number  of  segments.  As  we  will  see  later,  these  functions  are  especially 
expensive  to  approximate  using  uniform  segmentation.  The  fourth  column  labeled  Uni¬ 
form  will  be  discussed  in  Section  3.  The  fifth  (rightmost)  column,  labeled  Non-Uni/Uni, 
shows  the  ratio  of  the  number  of  non-uniform  segments  needed  compared  to  the  number 
of  uniform  segments  needed  as  a  percentage,  when  the  specified  approximation  error  £  is 
small.  For  example,  for  2X ,  non-uniform  segmentation  uses  84%  of  the  segments  needed 
by  uniform  segmentation.  For  y Fc  and  —  {x\og2x  +  (1  —  x)  log2(l  —  x)),  the  fifth  column 
contains  (0%),  which  shows  that  the  ratio  of  segments  in  non-uniform  segmentation  is 
vanishingly  small  compared  to  the  number  of  segments  needed  in  the  uniform  approxima¬ 
tion  as  £  — >  0.  We  discuss  the  derivations  for  the  number  of  segments  needed  for  uniform 
segmentation  in  Section  3. 


Table  1 

Number  of  Segments  for  Non-Uniform  and  Uniform  Segmentations 


Function 

f(x) 

Domain  of 

X 

Non-Uniform 

s  ~ 

Uniform 

s  ~ 

Non-Uni 
/Uni  = 

2X 

[0,1) 

0.2071 

W 

0.2451  _  In  2 

\/8Ve 

84% 

1/x 

[1,2) 

0.2071 

Ve 

0.3536  1 

y/s  y/Sy/e 

59% 

\fx 

[0,2) 

0.5946 

sTe 

0.0313  I 

e2  32s2 

(0%) 

l/^x 

[1,2) 

0.1378 

„  W 

0.2165  _  C3 
\fe  8  Vi 

64% 

log2(x) 

[1,2) 

0.2081 

W 

0.3003  _  I 

y/s  4-\/ln  2^/s 

69% 

111  X 

[1,2) 

0.1733 

Vi 

0.2500  _  I 

yfs  4  yfs 

69% 

sin(nx) 

[0,1) 

trims 

Ve 

U.o927  7 r 

Vi  8\/i 

76% 

cos  (it  x ) 

(0,1) 

072595 

vi 

0.3927  7r 

Vi  8\/i 

76% 

tan(  itx) 

[0,1) 

0.2005 

w 

0.3927  7 r 

Vi  8Vi 

51% 

y/-ln(x ) 

1  1  u 

0.5961 

6.9226  _  15fV/,(ln256)-2-i(ln256)-i 

9% 

L256’  4 2 

Ve 

Ve  V2  Vi 

tan2  {it  x)  +  1 

[0,1) 

0.4200 

yi 

0.7854  7 r 

Vi  _  4  Vi 

53% 

—  {x  log2  x+ 
(1-x)  log2(l  x)) 

(0,1) 

0.9434 

w 

0.2654  _  I 

e  2e In  2  e 

(0%) 

1 

[0,1) 

0.0550 

0.0754  _  Ve1+e_1  1 

73% 

l+e~x 

Vi 

Vi  e1+2+e— 1  4Vi 

l  -xZ 

^  2 

[0,V2] 

0.1452 

0.2233  _  V2 

Vi  4(27r)1/4Vi 

65% 
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3.  Uniform  Approximation  With  Unrestricted  Slope 

In  this  part,  we  consider  two  ways  to  determine  the  number  of  segments  needed  in  the 
case  when  a  completely  free  choice  of  slope  is  used  with  uniform  segmentation.  The  first 
approach  is  a  direct  computation,  which  can  be  applied  to  all  functions  considered  in  this 
paper,  and  the  second  is  an  asymptotic  approximation  that  applies  to  a  majority  of  the 
functions. 

The  first  approach  is  illustrated  as  follows.  Consider  a  uniform  segmentation  of  a 
function  f(x )  from  x  =  a  to  x  =  b.  Consider  a  segment  beginning  at  x  —  a  and  ending 
at  x  =  f3,  where  a  <  (3  near  a  point  in  the  domain  [a,b],  where  \  f'\x)\  is  maximum.  It 
is  at  this  point  that  the  maximum  error  between  the  function  and  its  linear  piecewise 
approximation  occurs.  That  is,  if  we  choose  the  segment  width  to  be  small  enough  at  this 
point  so  that  the  error  is  equal  to  the  specified  approximation  error  e,  then  that  small 
a  width  for  all  segments  will  be  sufficient  to  achieve  a  maximum  error  e  in  all  segments. 
Therefore,  the  number  of  segments  s,  to  achieve  an  approximation  error  with  uniform 
segmentation  is 


(8) 


Consider  a  piecewise  linear  approximation,  fp\(x)  to  f(x)  of  the  form  fp\(x)  =  ( f(/3 )  — 
/(a))|Ef  +  f(° 0-  The  error  due  to  the  approximation  can  be  viewed  as  |  f(x)  —  fp\\.  Note 
that  |  f(x)  —  /pi  |  is  0  at  x  =  a  and  at  x  —  (3.  However,  we  will  approximate  f(x)  in  the 
domain  [a,  /3\  by  adding  a  constant  to  fp\  so  that  the  maximum  error  in  the  domain  [ct,  /3\ 
is  no  greater  than  \\f(x)  —  fpi\-  As  a  result,  the  error  function  e(x)  for  the  domain  [a,  (3\ 
is  |(|/(x)  /pi | ) .  Substituting  for  fph  yields 


I  T  —  CM 

e{x)  =  x  lf(x)  -/(«)]  Hf(P)  -/(«)] -x — 

A  \_  p  —  a 


To  illustrate,  consider  the  function  f(x)  =  \fx  in  the  domain  [0,2],  f"(x)  becomes 
unbounded  near  x  =  0.  Thus,  we  choose  a  =  0,  and  from  (9),  we  have 


e(x)  =  ^  y/x-  • 

By  differentiating  e(x)  with  respect  to  x,  we  fold  that  the  maximum  error  occurs  at 
x  =  |.  At  this  value,  the  maximum  e,  emax  is  We  choose  this  value  to  be  e.  That  is, 
c max  =  £,  and,  so  ^  =  e.  Substituting  this  into  (8)  yields 


s  = 


1 

32^ 


(10) 
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Thus,  as  £  decreases  (improves),  the  number  of  segments  needed  for  a  uniform  seg¬ 
mentation  increases  as  the  inverse  of  the  square  of  e.  This  results  from  the  fact  that  f" 
becomes  unbounded  near  the  endpoint  0.  A  similar  analysis  can  be  applied  to  the  entropy 
function,  and  in  fact  can  be  generalized. 

The  number  of  segments  needed  for  uniform  segmentation  can  be  computed  for  all  of 
the  target  functions.  Column  4  of  Table  1,  labeled  Uniform  shows  the  expressions  of  the 
number  of  segments  required  in  a  uniform  approximation  for  the  various  functions.  For 
all  but  two  of  the  functions,  the  expressions  have  the  form  ^=,  which  are  the  same  form 
as  the  asymptotic  approximation  for  the  number  of  segments  required  in  a  non-uniform 
segmentation.  For  two  functions,  yfx  and  —  (x  log2  x  +  (1  —  x)  log2(l  —  x)),  the  asymptotic 
approximation  for  the  number  of  segments  is  -%  and  -,  respectively.  The  right  column  of 
Table  1,  labeled  Non-Uni/Uni  shows  the  ratio  of  segments  required  in  a  non-uniform 
segmentation  to  the  number  of  segments  required  in  a  uniform  segmentation  (expressed  as 
a  percentage).  For  some  functions,  like  2X,  sin( ttx),  and  cos(ttx),  this  is  high,  84%,  76%, 
and  76%,  respectively.  For  such  functions,  there  is  a  small  penalty  for  using  a  uniform 
segmentation.  For  other  functions,  like  yj —  ln(x),  this  percentage  is  low,  9%,  and  the 
penalty  is  high.  A  similar  statement  is  true  of  yfx  and  —{x\og2x  +  (1  —  x)  log2(l  —  x)), 
where  the  percentage  number  of  segments  is  expressed  as  (0%). 

ft  is  interesting  that,  when  the  domain  for  the  yfx  function  is  reduced  to  [1,2],  the 
number  of  segments  is  given  by  s  ~  °'^3. 

The  second  approach  to  determining  the  number  of  segments  requires  the  second  deriva¬ 
tive  f"(x )  to  be  bounded  over  the  domain  of  approximation.  We  have 

Theorem  2.  Consider  a  piecewise  linear  approximation  of  a  function  f(x)  on  the  domain 
[a,  b]  with  a  specified  approximation  error  e  or  less  using  uniform  segmentation.  Let  the 
absolute  value  of  the  second  derivative  \  f'\x)\  of  f(x)  on  the  domain  [ a,b ]  be  finite.  Then, 
the  number  of  segments  s  is 


(11) 


where 


(^-a)\Ar 


where  |/"|max  is  the  maximum  of  the  absolute  value  of  f"(x)  over  the  domain  [ a,b ]. 


Proof  See  Appendix. 

The  right  column  of  Table  1  shows  the  results  of  Theorem  2.  Specifically,  all  functions 
in  this  table  except  yfx  and  —  (log2x  +  (1  —  x)log2(l  —  x))  satisfy  the  restriction  that 
\f'(x)\  is  finite.  For  example,  for  sin( ttx),  \f"(x)\max  =  n2.  Therefore,  for  this  function, 
(12)  yields  s  ~  g^,  which  agrees  with  Table  1. 


4.  Concluding  Remarks 

As  a  result  of  our  analysis,  we  have  an  understanding  of  how  the  hardware  complexity, 
as  measured  by  the  number  of  segments,  depends  on  the  specified  precision  e.  Our  results 
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also  show  that,  for  some  functions,  it  is  reasonable  to  use  uniform  segmentation,  thus 
eliminating  the  segment  index  encoder.  For  such  functions,  our  results  validate  the  past 
research  on  uniform  segmentation.  For  example,  with  the  sin(nx )  and  cos( nx)  functions 
approximated  using  unrestricted  slope,  °"%27  segments  are  needed  for  uniform  segmenta¬ 
tion,  while  no  more  than  °'^95  segments  are  needed  for  non-uniform  segmentation.  This 
is  about  31%  more  segments.  The  penalty  is  substantial,  but  if  memory  is  inexpensive, 
and  speed  is  essential,  this  may  be  a  welcome  tradeoff. 

Table  2  shows  the  number  of  segments  needed  in  the  case  of  non-uniform  and  uniform 
segmentation  for  functions  approximated  using  unrestricted  slope  as  calculated  in  Sections 
4  and  5,  respectively.  Table  2  shows  the  number  of  segments  needed  for  four  precisions, 
8,  16,  32,  and  64  bits  3. 


Table  2 

Number  of  Segments  for  Non-Uniform  and  Uniform  Segmentation  For  Four  Precisions,  8, 
16,  24,  and  32  Bits. 


Function 

f(x) 

Inter¬ 
val  X 

8 

JNc 

16 

in-  U  nifc 
32 

)rm 

64 

8 

U 

16 

niform 

32 

64 

2X 

[0,U 

4 

75 

19,195 

1.26  x  10a 

6 

89 

22,717 

1.49  x  10a 

1/x 

[1,2) 

4 

75 

19,195 

1.26  x  10a 

8 

128 

32,773 

2.15  x  10a 

y/x 

[0,2) 

10 

216 

55,109 

3.61  x  10a 

8,206 

5.38  x  108 

2.31  x  1018 

4.26  x  108Y 

1/y/x 

[1,2) 

3 

50 

12,772 

8.37  x  10s 

5 

79 

20,066 

1.32  x  10a 

log2(U 

[1,2) 

4 

75 

19,228 

1.26  x  10a 

7 

109 

27,833 

1.82  x  10a 

111  X 

[1,2) 

3 

63 

16,062 

1.05  x  10a 

6 

91 

23,171 

1.52  x  10a 

sinfxx) 

[o,D 

5 

109 

27,759 

1.82  x  10a 

9 

143 

36,397 

2.39  x  10a 

cos(nx) 

(o,f) 

5 

109 

27,759 

1.82  x  10a 

9 

143 

36,397 

2.39  x  10a 

tan(nx) 

[0,1) 

4 

73 

18,583 

1.22  x  10a 

9 

143 

36,397 

2.39  x  10a 

y/-ln(x) 

y«m) 

10 

216 

55,248 

3.62  x  109 

157 

2,507 

641,600 

4.20  x  1010 

tali'1  (nx)  +  1 

[o,  U 

7 

153 

38,927 

2.55  x  10a 

18 

285 

72,793 

4.77  x  10a 

~{x  log2  x+ 

(1— *)  log2(l-a:)) 

(0,1) 

16 

342 

87,437 

5.73  x  10a 

136 

34,787 

2.28  x  10a 

9.79  x  1018 

1 

l+e~x 

[0,1) 

1 

20 

5,096 

3.34  x  10s 

2 

28 

6,989 

4.58  x  108 

—  ‘2. 

1  _ 

ysU  2 

[0,  y/2] 

3 

53 

13,458 

8.82  x  10s 

6 

81 

20,696 

1.36  x  109 

From  Table  2,  we  can  make  conclusions  about  the  feasibility  of  realizing  the  various 
functions.  Specifically,  we  can  see  that  for  64  bit  precision,  very  large  memory  size  is 
needed  in  all  cases.  For  32  bit  precision,  both  uniform  and  non-uniform  segmentation 
yield  feasible  realizations,  except  for  y/x  and  — (xlog2a;  +  (1  —  x)  log2(l  —  x)).  Also,  the 
memory  required  to  realize  a/—  ln(a;)  is  quite  large  compared  to  that  required  for  non- 
uniform  segmentation.  For  16  bit  precision,  all  realizations  of  the  functions  are  feasible, 
except  for  y/x  using  uniform  segmentation.  For  8  bit  precision,  all  realizations  are  feasible 
using  either  non-uniform  and  uniform  segmentation.  In  general,  for  many  functions, 

3 Assuming  that  the  most  significant  bit  is  the  coefficient  of  2  ,  we  choose  the  error,  e,  to  be  one-half 

of  the  value  of  the  least  significant  bit.  For  example,  for  8  bit  precision,  we  choose  the  error  to  be  2-9. 
We  substitute  this  for  £  in  the  equations  for  the  number  of  segments  for  non-uniform  segmentation  and 
uniform  segmentation. 
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uniform  segmentation  is  good,  especially  when  the  cost  of  memory  is  low. 
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5.  APPENDIX 

In  the  following  theorem,  we  make  the  following  reasonable  assumptions. 

Assumption  1:  As  £  — >  0,  the  measure  of  Ae  — >  0. 

Assumption  2:  The  open  set  B£  is  a  finite  union  of  open  intervals. 

Given  these  assumptions,  we  now  state 

Theorem  1  Consider  a  piecewise  linear  approximation  of  f  on  the  domain  [a,  b]  that  is 
accurate  to  within  e,  using  a  piecewise  linear  segmentation.  Let  f  be  three  times  con¬ 
tinuously  differentiable  on  [a,  6].  Then,  s(e),  the  number  of  segments  in  an  optimum 
segmentation  of[a,b],  satisfies  the  following  asymptotic  approximation: 

Q 

s(e)  ~  0), 


where 


c  =  \fa  V\f"(x)\dx ■ 

Proof  Given  e  >  0,  divide  the  domain  [a,  6]  into  segments  with  end  points  {xq,  x\,  . . . , 
where  xq  =  a  and  xs  =  b.  Assume  {xo,xi, . . .  ,2^}  has  the  fewest  segments  such  that  all 
segments  have  an  approximation  error  no  greater  than  e.  Thus,  for  any  segment,  if  we 
set 


Li{x)  =  f(x) 


f(xi+ 1)  ~  f(xi) 

Xi-\-i  Xi 


Xi)  +  f(x i)  , 


(13) 


then 


Li(x)  |  <  2e, 


for  Xi  <  x  <  Xi+i,  where  i  —  0, 1, . . . ,  s  —  1.  Without  loss  of  generality,  since  the  segmen¬ 
tation  is  optimal,  we  can  assume  that,  for  all  but  perhaps  one  segment,  there  exists  an  x* 
in  (x^  xi+i),  for  which  \L.fx*)]  =  2e,  for  i  =  0, 1, . . . ,  s  —  1.  Now  let 


( i )  T  Xi+ 1 

An  id  —  O 


and 
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Ni  Xi- {-i  Xi 

be  respectively,  the  midpoint  and  length  of  the  segment  [xi,Xi+ 1]  so  that 

(14) 

(15) 


This  yields,  after  some 

=  +  o(A»). 

Hence, 

M*)l  -  ^  I  A*l)l  (K2-1|  +  o(A,)), 

for  i  —  0, 1, . . .  s  —  1  and  Xi  <  x  <  Xj+i.  Since  —  1  <  a.i(x)  <  1, 


xi+ 1  =  + 


and,  for  x  G  [xj,a:j+i]; 


I  =  smid-Wy, 


where  —  1  <  <  1. 

Apply  Taylor’s  approximation  to  the  terms  in  L,  in  (13)  . 
algebra, 


Ae  =  {x  G  [a,  b]  :  \f"(x)\  <  \fe}  and 
B.  =  {x  £  [o,6]  :  |/"(l)|  >  -fi), 


(17) 

(18) 
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for  e  >  0. 

Thus,  we  split  the  sum  accordingly: 


s- 1 

i=  0 


A, 

V8 


A* 

y/8 


\/lf"(*(l)\  (1  +  O(A0) 

V7! f"(xl)\  (1  +  0(At)). 


The  first  sum  on  the  right  in  (19)  is  small,  0(£1/,‘t).  Also,  since 


(19) 


\f"(X^id)\  > 

for  xlmid  G  B£,  (16)  implies 

2c  >  A|/''(x«u)|(l  +  0(A,))  >  A^(1  +  o(A,)).  (20) 

It  follows  from  (20)  that  A)  =  0(y/e),  and  so  A,  — >•  0  as  e  — >  0.  Now,  (19)  can  be  written 
as 


s(e)  V2e 


0(«r1/1)+  ^ 

cl 


A* 


(1  +  O(Aj)) . 


(21) 


Since  A,  — >  0  as  e  — >  0,  (J  Z?e  =  [a,  6],  and  the  measure  of  A£  tends  to  0  as  e  — >  0,  it 
follows  that  the  number  of  terms  in  the  sum  in  (21)  goes  to  infinity  as  e  — >  0,  and  so,  by 
our  previous  assumptions,  we  can  approximate  the  sum  in  (21)  as  a  Riemann  integral: 


A) 

V8 


(l  +  0(Aj))  — 


1 

71 


(1  +  0(1)). 


Thus, 


s(e)  V2e  =  0(el/A)  +  ~^=  ^  yf\ f"(x)\dx  (l  +  o(l)) 

s(e)  =  0{e~1/4)  +  ~F~r  [  V\f"(x)\dx  (1  +  o(l)). 

j  Be 


and 
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Since  the  measure  of  the  set  A£  — >  0  as  e  — >•  0,  a  combination  of  the  above  results  gives 

Q 

s(e) - -=  ,,(e  -»•  0), 

V£ 

where 


yj\f"{x)\dx. 


Theorem  2  Consider  a  piecewise  linear  approximation  of  a  function  f(x)  on  the  domain 
[o,  b]  with  a  specified  approximation  error  e  or  less  using  uniform  segmentation.  Let  the 
absolute  value  of  the  second  derivative  \  f'(x)\  of  f(x)  on  the  domain  [a,  b]  be  bounded. 
Then,  the  number  of  segments  s  is 


s  ~ 


where 


(22) 


(b  -  "  M,  /" 


(23) 


where  l/'^max  is  the  maximum  of  the  absolute  value  of  f"(x)  over  the  domain  [ a,b ]. 
Proof  For  any  segment  in  a  uniform  segmentation  of  f(x),  the  difference  between  the 
exact  value  of  f  and  its  linear  piecewise  approximation  is 


Li(x)  =  f(x) 


f(Xj+ 1)  -  f(Xj) 
Xi- 1-1  Xi 


(x 


Xi)  +  f{Xi) 


We  require  that 
\Li{x)\  <  2e. 

By  a  process  similar  to  that  used  in  Theorem  1, 


(24) 


(25) 


L.(x)  =  \f"(x{2«)^ (a?W  -  1)  +  0(A?). 


(26) 


Let  /"|  assume  its  maximum  in  \xi,Xi+i]  at  x\.  Since  —  x*  +  O(Ai),  we  have 


l/"(4!d)l  =  l/>,*)l+0(Ai). 


(27) 


(25)  and  (26)  imply  that 
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max  |/"  0)| 

x€\Xi,Xi+i\ 


0a?M-l|+O(A?)<2£. 


(28)  can  be  written  as 


1 

lfe 


max  |/"  0)1 

x€[Xi,xi+i\ 


|a?(i)  -l|+0|0  <  A 


Now,  as  in  Theorem  1,  0  =  0(s  1//2),  so  that  (29)  can  be  written  as 

a9  1 

A2  < 


*  —  i 


ik  [maxxe[li,,i+l]  I/" 0)1]  010)  -  !|  +  O  (e"1/2) 


or 


(28) 


(29) 


A2  < 


16e 


[maxxe[xjjXj+1]  |/"0)|]  |a20)  -  !|  +  O  {e1/2) '  3° 

Now,  the  right  side  of  (30)  is  an  upper  bound  for  A2,  where  At,  i  =  0, 1, . . . ,  s  —  1  are 
s  segments  covering  [a,  b],  in  each  of  which  \Li(x)\  <  2e,  Xi  <  x  <  xi+\.  We  want  a 
uniform  segmentation  of  [a,  b]  with  the  minimum  Aj.  So,  we  choose  a,  =  0  m  (30), 
replace  [max^^^^^^  \f"(x)\]  by  the  maximum  of  \f"(x)\  over  the  entire  domain,  and 
take  the  square  root.  We  use  the  equality  sign  in  (30)  to  conclude  that 


A,= 


4\0 


^maXieH  LPOJK1  +  0(e1/2)) 


Thus , 


(31) 

=  l/"M  (6  -  a)  (1  +  0(e1/2)) 

Since 


s{e)  = 


b  —  a  dfe) 

A  i  sfe  ’ 


where 


lim  die)  =  d 
£— >0 


0 maxxe[aife]  I/" 0)1  (b  -  a) 
4 


we  see  that  (31)  implies  that,  for  a  uniform  optimal  segmentation 
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